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1.1 Introduction 

The use of spin polarized currents for manipulation of nano-magnetic structures, 
known as spintronics, has grown since it's first inception by Slonczewski and Berger 



( Berger, 1996 Slonczewski, 1996 ) into a thriving field of research ( Waintal et ai, 2000 



serkovnyak et al, 2002 Slonczewski and Sun, 2007||Ralph and Stiles, 2008[[Tserkovnyak et al, 2008[ 



Chudnovskiy et al., 2008||Bertotti, 2008|[Ralph and Stiles, 2009[ ). It has shown tremen- 
dous practical potential in it's ability to quickly read and write non- volatile informa- 
tion into high density memory structures by reversing the orientation of magnetic 
bits as well as it's ability to drive magnetic oscillators, which have a variety of uses 



dBurrowes et al., 2010[|Mizukami et ai, 2001[|Tatara and Kohno, 200"4l|Li and Zhang, 2004| 



Parkin et al., 2008 Tatara et al, 2008 Beach et ai, 2008 Hals et al., 2009 Tserkovnyak et al., 2005 



Silva and Rippard, 2008] [Cheng et al., 2010 ) . 

Slonczewski and Berger predicted that a current of spin polarized electrons passing 
through a ferromagnetic layer generates a torque on the magnet and thus can be used 
to change the orientation of the magnetization. This effect is known as a spin transfer 
torque. The spin-torque (ST) phenomenon consists of the transfer of angular momen- 
tum by the exchange interaction of electrons with the macroscopic magnetization. 
Due to the conservation of the total angular momentum in the interaction between 
spin polarized current and the magnetization of the ferromagnet, the net change in 
the angular momentum for the current flowing through the magnetic layer must be 
equal to the net angular momentum absorbed by the ferromagnetic layer. In turn, the 
angular momentum is proportional to the magnetization, the coefficient being given 
by the gyromagnetic ratio, hence the absorption of the angular momentum results in 
the change of the magnetization direction, that is in a torque on the ferromagnet. The 
existence of ST was confirmed experimentally ( [Myers et ai, 1999) [Tsoi et ai, 1998[ ) . 
This phenomenon provides the physical mechanism for the operation of spintronics 
devices. 

Magnetic nano-pilar structures typically consist of two magnetic layers, hereafter 
referred to as the "free" and the "fixed" layer, and one non- magnetic layer (e.g. a 
tunncling-transparcnt insulator or a metal) that are sandwiched in a pillar between 
two ohmic contacts. Fig. 11.11 The fixed and free layers differ significantly in their 
coercivity fields, with the fixed layer being pinned along it's orientation much more 
strongly than the free layer. This can be done by increasing the thickness of the fixed 



2 Dynamics of nano-magnetic oscillators 



layer, pinning it by proximity with an anti-ferromagnet, or picking a harder magnetic 
material. The free layer alternatively is much easier to manipulate. When an electric 
current is passed through the fixed layer in Fig. 11.11 in the ic-direction, it becomes 
polarized along the direction of the fixed layer. When it encounters the free layer the 
polarized spin-current induces ST and thus may change the magnetization direction 
of the layer, allowing for a number of dynamic regimes. 
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Fig. 1.1 Schematic diagram of magnetic tunnel junction (MTJ). Electrons run clockwise 
through the device first passing through the bottom "fixed" layer, where they become partially 
polarized in the — z-direction. They then pass through the "free" layer imparting some angular 
momentum. The MTJ has it's easy axis along the z-axis with the easy plane normal to the 
X-axis. 



One such regime involves using ST to switch the orientation of the free layer in a 
magnetic tunnel junction (MTJ). This use of ST was confirmed shortly after it's the- 
oretical prediction ( [Katine et ai, 2000] |Grollier et ai, 2001[ ) with much of the early 



research focusing on understanding the effects of internal anisotropies (Sun, 2000 
|Mangin et ai, 2006[ ) , temperature ( Ozyilmaz et al., 2003| Urazhdin et ai, 2003, Hosomi et ai, 2005[ ) , 
and spin-current strength ( [Myers et al, 2002 Seki et al., 2006 1 on the switching dy- 



namics. With the development of high tunneling magnetoresistance MTJ's ( Yuasa and Djayaprawira, 2007) 
and with switching times in the sub-nanosecond regime dTulapurkar et ai, 2004||Kent et al, 2004D 
ST random access memory has become a very attractive candidate for non-volatile 
memory. While development of ST memory has moved forward ( [Matsunaga et al., 2009P 
recent research has also focused on more exotic switching strategies. Some of these 
strategies include using high frequency ac spin-currents to excite resonant spin modes 
to assist the switching (Cui et ai, 2008) as well as using spin-current polarized perpen- 
dicular to the interface ( [Nikonov et al., 2010[ ) to increase energy efficiency and reduce 
switching times ( [Bedau et ai, 2010 ). 

Another major area of research for MTJs focusses on the use of MTJ's as small 
oscillators, known as spin-torque oscillators. Oscillations in MTJ's were first observed 
in some of the earliest studies of ST (Kiselev et ai, 2003 Rippard et al., 2004). Since 
ST oscillators can operate in GHz frequencies they've been an attractive candidate 
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for use in timing mechanisms for computers and in radio-frequency devices for use 
in radar and telecommunication devices ( |Deac et ai, 2008] ). In the last few years 
it has also been shown that ST oscillators can exhibit resonant excitations with 



ac spin currents which oscillate near their natural frequency (Georges et ai, 2009 



Ruotolo et ai, 2009| pieng et ai, 20W\ [Finocchio et ai, 2010| [Rippard et ai, 2010 



Finocchio et ai, 2011 1. Because of these practical applications, much of the research 



done in the last decade has centered on the reduction and control of the oscillation 



linewidth (Mizukami et ai, 2001 Krivorotov et al, 2005 Nazarov et ai, 2006 Pribiag et ai, 20091 



While it's been shown that external fields (Thadani et al., 2008 Braganca et al., 2010 
and oscillation amplitude ( |Rippard et al, 2006 Kim et al, 2008| Boone et al, 2008 



can be used to tune the linewidth, the linewidth itself is the result of noise in the 
MTJ and much research has looked at understanding sources of noise and decay 
in ST oscillators ( |Sankey et al., 2005| [Mistral et al., 2006| [Tiberkevich et al, 2007| . 
Thus understanding of stochastic processes such as thermal and shot noise has be- 
come increasingly important for operation of ST devices, as manufacturing technology 
has improved and these devices have been pushed to smaller and smaller sizes. 

First description of thermal noise in ferromagnets has been developed in the pio- 
neering work by W. F. Brown ( [Brown, 1963[ ), where the noise have been described as a 
stochastic magnetic field acting on the magnetization. The Fokker-Planck equation for 
the probability distribution of magnetization direction has been formulated, and served 
as a basis for consideration of the noise in various dynamical regimes. Further devel- 
opment of theoretical treatment of the noise has been provided by IForos et al. 2005] 
IDuine et al. 20071 IForos et al. 20081 A crucial simplification of the theoretical descrip- 
tion of the magnetization noise has been reached in the work of [Apalkov and Visscher 2005 



by noticing that the dynamics of ST device exhibits a time scale separation between 
the fast magnetization precession and slow change in the energy of the precessional 
orbit. This allowed to average the equations of motion over the fast variables and 
formulate the one dimensional Fokker-Planck equation for the energy distribution. 

In this chapter we explore how stochastic processes effect switching elements as 
well as ST oscillators. To do so we first discuss the deterministic dynamics of MTJs 
followed by the corresponding stochastic dynamics. Our particular emphasis is on the 
effects of out-of-equilibrium noise, such as spin-torque shot-noise, inherent to any de- 
vice operating under an applied current. Following Ref. ( [Apalkov and Visscher, 2005[ ), 
we derive a stochastic Langevin equation for the slowly varying energy of precessional 
orbit, thereby introducing the energy noise and the energy diffusion coefficient. We 
then formulate a Fokker-Planck equation for the energy density distribution. Wc use 
these constructions to analyze switching time distribution as well as the shape of the 
optimal spin-current pulse, which minimizes Joule losses of a switch. Wc also derive 
generic expression for the linewidth of ST oscillator and discuss its dependence on 
temperature, spin-current amplitude and other parameters. 



1.2 Deterministic dynamics 

In this section we briefly discuss magnetization dynamics of a mono-domain ferromag- 
net with the emphasis on the time scale separation, which allows to introduce magnetic 
energy as a slow (hydrodynamic) variable. On the most basic level the magnetization 
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M precesses around the instantaneous effective magnetic field Hcff according to the 
Landau-Lifshitz (LL) equation 

MLL = -7[MxHeff], (1.1) 

where 7 is the gyromagnetic ratio. The effective magnetic field is given by the gra- 
dient of the magnetic energy i?(M) with respect to the magnetization, i.e. Hcs ~ 
— Vm^'(M). The energy of magnetic layers, conventionally used in magnetic junction 
devices, typically includes a combination of an easy plane and the in-plane easy axis 
anisotropy as well as the effect of an externally applied magnetic field 



E{M) = fio 



-g-(M.zf + ||(M..f-H...M 



(1.2) 



Here represents the strength of the easy axis anisotropy, chosen in the z-direction, 
while represents the strength of the easy (y— z)-plane anisotropy field. Fig. II. II The 
anisotropy fields create an energy landscape with two valleys separated by an energy 
barrier. The minima of these valleys determine the orientation of the two stable states. 
External fields Hoxt arc often used to force the magnetization to align in a specific 
direction or to allow one of the energy minima to be lower than the other. 

The Landau-Lifshitz equation (|l.ip evidently conserves both the absolute value of 
the total magnetization |M| = Ms (subscript s stays for the saturation magnetization) 
and the magnetic energy i?(]V[). As a result, the magnetization is confined to move 
along a closed orbit of a constant energy £'(M) — E, belonging to a sphere of radius 
Ms- These are the so-called Stoner-Wohlfarth (SW) orbits, plotted in Fig. 11.21 




Fig. 1.2 Several closed SW orbits of constant energy in the northern and southern hemi- 
spheres for HI = HI = Ms. 

It was realized long ago that the actual magnetization dynamics must include the 
energy dissipation, originating from the coupling of the global magnetization M with 
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the host of microscopic degrees of freedom. The simplest phenomenological way to 
incorporate these effects in the magnetization dynamics was suggested by Gilbert. 
The corresponding Gilbert damping (GD) term in the equation of motion leads to the 
decay of the precession and aligns the magnetization along the effective magnetic field 

Mgd = ^[MxM]. (1.3) 

The strength of the Gilbert dissipation is proportional to the phenomenological di- 
mensionless damping constant a. In modern nano-magnetic devices its value may be 
as smah as a = 0.01 ( |Myers et al, 1999[ [Katine et al, 2000[ [Ozyilmaz et al, 2003D , 
allowing for dozens precession cycles prior to equilibration. 

Dramatic progress of the last decade is associated with the realization by Slon- 
czewski ( Slonczewski, 1996) and Berger ( |Berger, 1996[ ) that the magnetization may 



be influenced by the spin polarized current of electrons through the ST effect. Enter- 
ing the free layer, spin polarized electrons found themselves either aligned (with the 
amplitude oc cos 0/2, where 9 is an angle between polarizations of the fixed and free 
layers), or antialigned (with the amplitude cx sin 0/2) with the free layer magnetization 
direction. In the latter case the strong exchange interactions cause the rapid spin-flip, 
transferring the angular momentum h from the itinerant electron to the macroscopic 
polarization of the free layer. The corresponding non-conservative term in the macro- 
scopic equation of motion takes the form 



M X [is X M] , (1.4) 



which tends to align the magnetization along the spin polarization of the electric 
current. The spin-current vector Ig is directed along spin polarization of the fixed 
layer, Fig. ll.ll Its absolute value is proportional to the difference of the electric currents 
of spin-up and spin-down electrons, Xg = (J-j- — I^) fi/2e. 

The three above mentioned terms provide a conventional deterministic description 
of the magnetization dynamics 

M = Mll + Mgd + MsT. (1.5) 

Fig. 11.31 shows magnetization vector as well as the azimuthal angle 6 and magnetic 
energy (|1.2p as functions of time in the regime corresponding to the current-induced 
magnetization switching experiment. In this case the spin current is sufficiently strong 
to overcome the damping and make magnetization to cross the energy barrier, eventu- 
ally coming to rest along negative z-axis. As Fig. ll.Bl shows the magnetization precesses 
around the z-axis many times before crossing over the energy barrier. It may be no- 
ticed that, while the angle oscillates widely on the time scale of the precession period, 
the energy is a smooth function, which changes only at a much longer time scale. The 
reason for this behavior is that the energy-conserving LL term in eqn (|1.5p is typically 
much larger than energy-non-conserving dissipative GD as well as driving ST terms. 

This observation suggests to parameterize the instantaneous magnetization direc- 
tion by the slow energy variable E and a fast 27r-periodic variable which runs uni- 
formly along a closed SW orbit of a constant energy. Up to precession period Pe the 
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Fig. 1.3 (Left) Switching trajectory for = H% = Ms and Is = -0.05A/s. (Right) Az- 
imuthal angle 6 vs. time (full line) and energy vs. time (dashed line) for magnetization reversal 
of the system with a = 0.01, = 0.028Ms, = M^, and Is = -O.OlAf^. Time is measured 
in units of to = (^Ms)^^ with energy in units of tvEq = 

latter is simply the time needed to reach a given point along SW orbit dif = 27r dt/Pg;, 
where according to Eq. (jl.ip 



dt 



dM ■ [Heff X M] 
7|[HeffXM]|' 



Pe 



dM ■ [Hoff X M] 

7 I [Heff X M]|' 



(1.6) 



Here the integral runs along the orbit of energy E. We thus notice that the parametriza- 
tion M = M(i?, (p) is specified implicitly by the following relations: 



1 



[Heff X M] 



9bM 



[M X [Heff X M]] 
I [Heff X M]|' 



(1.7) 



where VLe = Ztt/Pe is the precession frequency. The first one is a consequence of LL 
equation (|l.ip and our choice of while the second can be checked by e.g. going 
to the instantaneous reference frame, where Heff is directed along the z-axis. Notice 
that d^pM ■ Q^M = and thus the new coordinates are locally orthogonal. Such 
parameterization can't be global, of course, but rather should be introduced separately 
for each of the four basins apparent in Fig. 11.21 Notice also that 



dAE 
dE 



r M ■ [OeM X dM] / [M X OeM] ■ dM 
f K K =^^^^^-^ 



(1.8) 



where Ae is the oriented area of the sphere, enclosed by the orbit with energy E. Here 
we employed eqs (|1.7p and (|1.6p . Equation (|1.8p provides geometrical interpretation of 
the period as a change of the orbit's area. 

One can write now the equation of motion (|1.5p as equations for i?(t) and (pit). To 
this end we notice that 



E = Vm-E • M = -Heff • M ; 



Qe [Heff X M] ■ M 
7 |[HeffxM]|' ■ 



(1.9) 
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and substitute M from eqn (jl.Sp . The fast LL part obviously drops from the equation 
for the energy, while the remaining slow terms acquire the form 

E ^ -F{E,^)+i,-ViE,if). (1.10) 

The two generalized forces on the right hand side originate from the Gilbert damping 
and spin-torque 



a 



1 



F{E,ip) = — MxHeff' -M; V(i;,^) = ^ MxM . (1.11) 



correspondingly. Hereafter M is understood as a conservative LL part, eqn (|l.ip . only. 
On the other hand, Gilbert damping drops from the equation for the angular variable 
Lp and its dynamics is mostly governed by the uniform precession 

^^n. + n /^-^^'^^i . (1.12) 

The second term here originates from the fact that the spin-torque may have a com- 
ponent along the SW orbit. This leads to a certain renormalization of the precession 
frequency Qe- Since the integral along the closed orbit of the second term on the right 
hand side of eqn ()1.12|) vanishes, such renormalization is of the order 0{i^). In the 
spirit of time scale separation wc shall neglect this kind of corrections. Similarly we 
disregard all terms which arc of the order 0(a^), assuming a small Gilbert damping 
constant. 

The above remarks show that the fast variable Lp approximately exhibits the uni- 
form precession according to 93 = VIe = 2tt^ MgdE / dA, cf. eqn (|1.8p . allowing to 
identify A/2TT'-fMs with the action canonically conjugated to the angle (p. It suggests 
to average equation (|1.10[) for the slow energy variable over the precession period, 
arriving at 

E = -FE+is-yE (1.13) 
where the averages of the generalized forces are given by 

FE = Y7ir {'^[dMxU,s]-M; = -rr^ / [rfM x M] (1.14) 

MsfE J M^fE J 

correspondingly. The integrals here run along SW orbit with energy E. Notice that to 
evaluate these forces one does not need any information about dynamics, but rather 
only the form of static SW orbits. While this strategy was shown useful in the analysis 
of non-linear oscillators long ago ( [Dykman and Krivoglaz, 1979 1 it was probably first 



applied in the present context by Apalkov and Visscher ( Apalkov and Visscher, 2005| 
and further developed by Bertotti ( |Bertotti, 2008[ ). 

Figure fTT^ shows the total generalized force as a function of energy for various values 
of the spin-current I, • For sufficiently small spin-current the force is strictly negative 
(dashed line, E < Q) forcing the magnetization to relax to the equilibrium energy 
minimum. At a larger current there is an energy value i?, where E = Q, correspond- 
ing to a stable precession around SW orbit with the energy E ( [Kiselev et al., 2003{ 
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[Mistral et al, 2006D . Finally at a spin current larger than a certain critical value the 
generalized force is positive in the entire energy interval (dot-dashed line, i? > 0). As a 
result the magnetization is forced to switch to the new stable position. In the following 
sections we examine effects of noise on both switching and stable precession. 




0.0 0.2 0.4 0.6 0.8 1.0 
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Fig. 1.4 Shows the generalized force, E = —Fe + is ■ Vb vs. energy for various strengths 
of Is for system with = 0.028Ma and = Ms. From top to bottom Is = O.OOSAfs 
(dot-dashed), Is = 0.006Afs (solid). Is = 0.004J\fs (dashed). 



1.3 Equilibrium stochastic dynamics of magnetization 

The deterministic description of the previous section is incomplete, especially for 
small enough magnetic domains. The reason is the stochastic part of the magnetic 
torque. The fact that noise must accompany Gilbert damping to satisfy the equilibrium 
fluctuation-dissipation theorem (FDT) was first realized by W. F. Brown ( [Brown, 1963[ ) 
The way he introduced the noise was by adding a stochastic component to the effective 
magnetic field in eqn (jl.ip . Hcff — > Hcs + h(i), where h{t) is an isotropic Gaussian 
noise. Its amplitude may be uniquely determined from FDT and is given by 

{h'{t)h^{t')) = 25^j5{t - t')akBT/M,-f = 2^^J6{t - t')D, (1.15) 

where T is the temperature and D is diffusion constant. In presence of the noise, 
the magnetization dynamics should be described by a probability distribution of the 
magnetization vector ^(M,^), obeying a proper Fokker-Planck (FP) equation. 

We shall first outline how it works in an equilibrium setting, generalizing it to the 
presence of the non-equilibrium spin-current in the next section. Stochastic magnetic 
field adds a Langevin force term Mstoch = 7 [h(f) x M] to the equation of motion ()1.5|) 
for the magnetization. It leads to Langevin force terms in equations of motion (|1.10|) . 
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(|1.12p for the energy-angle variables. The explicit form of the Langevin forces can be 
obtained using the projection on ip- and i?-directions respectively, as it was done for 
deterministic eqs (|1.9p . As a result wc obtain the stochastic equations of motions with 
the multiplicative noise 

E = -F{E,ip)+SE{E,^)-h{t); (1.16) 
ip = nE + g4E,ip)-h{t), (1.17) 

where the two mutually orthogonal noise-multiplying vectors are given by 

gE^M^nEd^M; =7nB[M X M]/|M|2 = -riB^BM (1.18) 



and we employed eqs (|1.7p . Since the multiplicative noise originates from the change 
of variables M E,ip, it should be treated with the Stratonovich regularization 
( [van Kampen, 2dOT| . Employing the orthogonality • gi^ =0, one obtains the fol- 
lowing FP equation for the probability density V{E, ip, t): 



V^df 



{F-kE)V + DdE[\gE\'V) {-VLe-K)V + Dd^[\g^\'V 



|2 , 



(1.19) 

where the diffusion coefficient D = ak-^T /Mg^ is determined by the equilibrium cor- 
relator of stochastic magnetic field eqn (|1.15p . The "false" forces fc^ and k^p are given 
by fcp = D^u{dv^pi) ■ ( |van Kampen, 20dT] |Kamenev, 201 1[ ), where = E,(f. 
Employing eqs (|1.18p along with orthogonality 9<^M • Q^M = 0, one finds 

kE = ^dE (nl |9^M|') + ^nldE \d^M\^ = D{dE - dEln^E) |g£|' ; (1.20) 

k^ ^ (n% laeMi") + ^iild^ la^Mi' - Dd^ |g^|' . (1.21) 

This transforms the FP equation (|1.19p into 

V^dE [{F + VeOe \nnE)V + VeOeV] + [-VLeV + V^d^V] , (1.22) 

where the two diffusion coefficients are 'D^{E, ip) = D |gp|^. 

To proceed we employ the periodicity of V{E,ip,t) in the angular directions to 
write it as a Fourier series V{E,ip,t) = X]m=-oo ^ni(^' 0^'™"^- ^"^^ expects that on 
the fast time scale of the precession period Pe the distribution nearly equilibrates in 
the angular direction and then evolves slowly along the energy direction. This implies 
that at long times t ^ Pe one has Vq ^ Vm, where m 7^ 0. As a result one can write 
a closed FP equation for Vo{E,t) in the following form 

= dE [{Fe + VeOe \nnE)ra + VeOeVo] , (1.23) 

where the average generalized force Fe is given by cqn (|1.14p and the averaged over 
the period energy diffusion coefficient is 

'De = ^ DVt\ |9>pM|^ dip = ^ j>DM. ■ dM ^ ^ j>D[dM. y. Hoff] • M. (1.24) 

One may now include V±i, V±2, etc. as progressively small corrections and take into 
account how the incomplete angular equilibration affects evolution in the energy di- 
rection. We shall not take this root here. 
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Notice that in equilibrium FDT dictates D = ak-gT/Msj and thus T>e = k-QTFE, 
cf. egn (jl.l4p . As a result of this remarkable relation the equilibrium stationary solution 
of the ID FP equation (|1.23p takes the universal form 



Vq[E) = i exp 



Pe 
2ttZ 



exp 



E 
k^ 



(1.25) 



where Z is a normalization factor. The very last term here is, of course, the Boltzmann 
exponent. It is crucial that the entropy S{E) = fceTlnrig^ has emerged, which may 
be traced back to the density of states along SW orbit given by Pe oc dA/dE, cf eqn 
(|1.8p . Notice that it is the "false" force kE, eqn (|1.20p . which is responsible for this 
fact. 



1.4 Shot-noise and non-equilibrium stochastic dynamics 

The advent of spin-polarized currents put the operation of spintronics devices in far 
from equilibrium conditions. This not only influences the deterministic motion through 
the ST term ()1.4|) . but also changes the stochastic force. Unlike the equilibrium noise, 
where the noise is uniquely determined by FDT, the out-of-equilibrium noise is non- 
universal and depends on details of a specific setup. Here we consider nonequilibrium 
noise in a spin valve device consisting of two magnetic layers separated by a tunneling 
barrier, the so-called magnetic tunnel junction (MTJ). The MTJ is modeled by the 
two itinerant ferromagnets, whose majority {a = +) and minority {a = —) bands are 
described by the operators cj,^,Cka for the fixed ferromagnet and dj^,dia for the free 
layer. The corresponding Hamiltonian consists of three parts describing isolated fixed 
and free layers along with the tunneling of electrons between them 

-f^ = -f^fixcd + ^^froc + -fftun- (1-26) 

The fixed layer is described by a Fermi-liquid Hamiltonian 

^fixod = ^ ekacl^Cka- (1-27) 

The local exchange interaction of itinerant electrons with magnetic momentd in the free 
layer can be written as Jq J drS{r) ■ s(r), where S(r) = M(r)/7 is the spin density of 
the magnetic layer and s(r) denotes the spin density of itinerant electrons. We assume 
that the spin current flowing through the free layer does not disturbe its monodomain 
magnetic structure. In that case, the layer can be considered as a single large spin 
with the magnetic moment Mg. The Hamiltonian of the free layer that accounts for 
the interactions of the itinerant electrons in the free layer with its total spin S ~ Ms/7 
is given by 

Htree = J2 - JS • S - 7S • Heff , (1.28) 

la 

where s = \ Yl,iaa' ^l^^aa'dia' is the spin of itinerant electrons and J is Heisenberg 
exchange interaction constant. Finally, the tunneling term in the Hamiltonian is given 

by 
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Htun = w^r'^lA'^' + (1-29) 



Here the spin indices of the operators and dia' denote the spin projections along 
the magnetization directions in the fixed and in free layers respectively. Because of 
a finite angle between the two magnetizations the tunneling matrix elements become 
spin-dependent. They are given by Wl^f = {(7\a')W, where the spin-transformation 
matrix is {a\a) = e"*'"'^/^ cos 6*72 and {a\a') = e*'""^/^ sin 6*72. The angles {0, (f) are the 
polar and azimuthal angles that denote the direction of the magnetization in the free 
layer in the reference frame with z axis pointing in the direction of magnetization of 
the fixed layer. 

The main source of the non-equilibrium noise is the discrete nature of the angular 
momentum transfer. Indeed, each electron tunneling at a random time may undergo 
a spin-flip with a probability depending on the mutual orientation of the quantization 
axis. Since each event transfers exactly the unit of angular momentum /i, the direc- 
tion of an ensuing magnetization rotation is random due to the uncertainty principle 
(i.e. non-commutativity of the three components of the angular momentum operator). 
Therefore the electron transport is accompanied by the stochastic force acting on the 
magnetization. Its nature is similar to the charge shot noise. While the latter is due to 
the discreteness of charge, the former is due to the quantization of the angular momen- 
tum. We call it thus the spin shot noise ( Foros et ai, 2005] [Chudnovskiy et ai, 2008 



Swiebodzinski et ai, 2010) 



To introduce this noise into the semiclassical equation of motion we write the model 
specified by eqn ()1.26|) as a path integral on the Keldysh contour ( [Chudnovskiy et ai, 2008} 
[Swiebodzinski et ai, 2010} [Kamenev, 201 1[ ) and integrate out all fermionic degrees of 
freedom, while keeping the dynamics of the macroscopic spin S. The integration is 
facilitated by the tunneling approximation, i.e. keeping only the lowest nonvanishing 
terms in the coupling W. As a result the spin torque and spin shot noise terms as well 
as an additive renormalization of the Gilbert damping are generated. Furthermore, the 
spin shot noise can be cast into the form of random fluctuations of the spin current 
vector 5Is entering the spin torque term p.4p . The noise correlator can be expressed 
in terms of MTJ conductance for parallel Gp and antiparallel Gap magnetization 
orientation of the two layers, (Gp > Gap) 



where 



(SmSii(t'))=2M;D(l))S,j5(t-t'), (1.30) 

where ao is a bare Gilbert damping of an isolated layer and V is a voltage bias between 
the two ferromagnets. The first term in eqn (|1.3ip describes the equilibrium noise due 
to intrinsic relaxation processes in the active layer, it is therefore proportional to the 
intrinsic Gilbert damping constant ag. The second term describes the additional noise 
due to the coupling of the active layer to the fixed layer which in this case plays 
the role of the external reservoir of spins. The external contribution to the noise is 
proportional to the spin- flip current Isf (6*). The latter counts the total number of spin 
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flips irrespective of the direction of the ensuing magnetization change (as opposed to 
the spin current I5). In the magnetic tunnel junction setup we found for the spin- flip 
conductance 

disfie) 

— I T D ^^l r I I 

2, 



dV 

Gp = G-i 



h 
4^ 



Gap cos^ 



(1.32) 



G — ; Gap — G-]-— + G— 



+ ) 



where we adopted notations of ( [Slonczewski and Sun, 2007] ) for the partial conduc- 
tances Ga-a' between the spin-polarized bands of the two ferromagnets. The external 
spin noise is accompanied by the renormalization of the Gilbert damping constant, 
which now depends on the instantaneous magnetization of the free layer through the 
angle 9 it forms with the magnetization of the fixed layer 



a{9) = ao + 



h-f f duje) 

dV 



= a(M) , 



(1.33) 



where the spin- flip differential conductance is given by eqn (|1.32p . For low voltages, 
eV < k^T, the tunneling contribution to the noise is proportional to the temperature 
and the total noise satisfies FDT with the renormalized Gilbert damping constant 
a(M), ( |Foros et ai, 2005] [Chudnovskiy et ai, 2008D . For voltages exceeding the tem- 
perature eV > k^T, the external noise is essentially nonequilibrium, proportional to 
the applied voltage and independent of temperature. 

The conductances Gaa', Gp, and Gap allow the complete characterization of the 
electric and magnetic properties of MTJ. So, the electric conductance of the MTJ in 
an arbitrary orientation is given by dig {6)/ dV ~ Gp cos^ {9/ 2) + Gap sin^ {9/2). Notice 
that the spin shot noise is minimal for P orientation and maximal for AP one - exactly 
opposite to the charge current and the charge shot noise. In contrast to the spin-flip 
current introduced above, the spin-current Xg is governed by the spin conductance 



(Slonczewski and Sun, 20071 



dXs 
'dV 



4e 



(G^ 



G— 



G+- -G__). 



(1.34) 



Taking noise into account, one arrives at the stochastic equation of motion for 
magnetization 



M = ~7 [M X Heff] 



«(0) 
Ms 



M X M 



7 
M2 



M X 



is + sut)) X M 



(1.35) 



The equivalence of the two forms of description of the magnetization noise, cither as 
a stochastic spin current or as a stochastic magnetic field, is expressed by the relation 
of the corresponding cumulants 



^ {Smsiiit')) ^ 2D{M) % S{t - t') . 



(1.36) 



Notice that the noise correlator itself (|1.3ip as well as the renormalized Gilbert damp- 
ing (|1.33|) depend on the instantaneous angle between the free and fixed layers magne- 
tization direction. Being derived from the path-integral representation, the correlator 
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(jl.36|l should be regularized in the retarded (Ito) way. Moreover, in presence of the 
spin-torque the Langevin equation ()1.16|) for the energy variable acquires the addi- 
tional deterministic term Ig • 'V{E,ip), see eqs (|1.13p and (|1.14p . This term describes 
pumping of the energy into the magnetic system by the spin torque process. Together 
with the Gilbert damping, the spin torque term determines the dynamics of the spin 
switching process as well as the average energy of the SW orbit for the steady state 
magnetization precession. 

One can now generalize derivation of FP equation (jl.23|) . presented in section [TT3l 
for the presence of the spin-torque and non-equilibrium noise. The result is given by 
the continuity relation in the energy direction Vo + Be Jo = 0, where the probability 
current in the energy direction is given by 



ME,t) = - 



- Is • V, 



VeOe Infii 



Vq - VeOeVo ■ 



(1.37) 



Here the components of the deterministic force and the energy diffusion coefficient are 
given by eqs (|1.14p and p.24p . correspondingly, where one has to substitute a — > q;(M) 
and D — > D(M) in accordance with eqs (|1.33p and (|1.3ip . 

The stationary solution of FP equation (|1.37p plays a special role in several regimes, 



such as activated magnetization switching ( Apalkov and Visscher, 2005 1 or steady 
state magnetization precession ( Kiselev et ai, 2003 ), considered below in more details. 
This solution is obtained by putting j7o = in eqn (|1.37p and reads as 



VoiE) cx Pficxp 



Fe -Is -Vi 
Ve 



dE 



(1.38) 



Notice that in equilibrium, i.e. Ig = and D(M.) = a(M)fcB'r/Ms7, one has "De = 
k^TFE, even for the magnetization-dependent damping constant and the noise cor- 
relator! As a result the Boltzmann form (|1.25p with the entropic factor given by Pe 
still holds. However, away from equilibrium the stationary energy distribution (jl.38p 
may be rather different from the Boltzmann shape. The ID Fokker-Planck equation 
(|1.37p is especially useful to analyze non-stationary regimes of spin- valves, such as e.g. 
kinetics of the switching process. 

1.5 Dynamics and optimization of the magnetization switching 

According to the deterministic equation of motion (jl.l3p one may introduce the energy- 
dependent critical spin-current 

Fe 



IciE) = - 



(1.39) 



which nullifies the right hand side of eqn (|1.13p . The critical spin-current is then 
defined as Ic = ma.x{Ic{E) | < _E < Eq}, so if the applied spin-current Is > Ic it 
forces the magnetization direction to switch to a new stable energy minimum. The 
corresponding deterministic switching time is given by 



(Eini) 



dE 



e,^,z-Ve{Is-Ic(.E)) 



(1.40) 
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Fig. 1.5 Numerical simulations of the switching process are shown for — 0.028Afs, 
= Ms, and T = SOOK. (Left) Switching probability, as a function of time for spin-current 
values 1.61c (dot-dashed), l.llc (solid), and 0.8Ic (dashed). (Right) Average switching time 
as a function of spin-current relative to the the critical current, Xc. Time is measured in units 

of (7Ms)"' 

where is an initial energy. This expression diverges as (Is — Ic)^^ and such a 
tendency was indeed observed in experiment ( |Bedau et ai, 20 IQ] ). This divergence is 
augmented by the stochasticity, since even for the currents somewhat less than critical 
the switching docs occur, albeit taking exponentially long waiting time. The proper 
description of the switching time must therefore rely on the probability distribution 
to undergo the irreversible switch during time tgw Such probability distributions cal- 
culated by Monte Carlo simulations of cqn (|1.35|) arc presented in Fig. 11.5b . for several 
values of the spin-current Xg. Having such distributions, one may evaluate e.g. the 
average switching time, plotted in Fig. 11.5b as a function of the applied spin-current. 
One observes that the switching time indeed grows exponentially at Is < Ic- An- 
alytically this latter regime was analyzed using eqn ()1.37p by Apalkov and Vissher 
( [Apalkov and Visscher, 2005t . Although Ref. ( Apalkov and Visscher, 2005t used the 



equilibrium noise, we found ( [Swiebodzinski et al, 2010 ) that inclusion of the shot 



noise (|1.31l) merely changes the effective temperature and therefore is hardly distin- 
guishable from effects of heating. 

One may notice that eqn (|1.40|) for the deterministic switching time also tends to di- 
verge if -Eini "> 0. This is due to the fact that if the initial direction of the magnetization 
is very close to the magnetic energy minimum, it formally takes a logarithmically long 
time to depart from the minimum. Clearly such a spurious divergency disappears once 
the distribution function of the initial (i.e. prior to the application of the spin-current) 
magnetic energy is taken into account. It is reasonable to assume that such a distribu- 
tion is thermal and is given by eqn ()1.25|) 7'o(£'ini) oc Pei^i exp [—Eini/kBT]. Although 
the (non-equilibrium) noise acts on the magnetization during the entire switching tra- 
jectory, we found that the simulation data for > Ic, Fig. 11.51 can be rather well 
fit by taking into account this initial distribution only; i.e. assuming the deterministic 
motion with the random initial condition p(isw) = / cJ^-ini 7'o(£'ini)<5(tsw ^ ^sw(-E'ini))- 
We therefore conclude that measuring the switching time distribution. Fig. 11.51 is not 
a sensitive way to resolve the effects of the non-equilibrium spin shot noise. 
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The treatment of the stochastic dynamics based on the time scales separation, pre- 
sented in the previous section, is nevertheless very useful to analyze an optimal pro- 
tocol for the time-dependent spin-current pulse Is{t). The issue at hand is as follows: 
as shown above, to make the switching time shorter one needs a larger spin-current. 
This leads to an excessive Joule losses cx / dtl^{t). On the other hand, decreasing the 
spin-current leads to a longer switching times, which in turn may again increase heat- 
ing due to the long time needed to complete the switch. One may expect that there 
is an optimal strategy, which minimizes the heating losses. To find such a strategy 
we treat the system as a heat engine and require that the ratio of the instantaneous 
energy gain E(t) to the energy loss oc Ig{t) is maximized at each instance of time. 
Substituting the energy equation of motion (|1.13|) and taking variation with respect 
to the spin-current, wc find the following condition 

Solving for the spin current , yields 

I°P* = 2,,^ = 2Ie(^;), (1.42) 
z • V_B 

cf. eqn (|1.39p . Therefore the optimal current is exactly twice the one needed to nullify 
the energy flow in eqn p.l3|) . Substituting this back into eqn (|1.13p . one finds that 
upon the optimal current the latter takes the form 

Fe. (1.43) 

This is exactly the same as in eqn (jl.lOp without the external spin current Is, but 
with the time being reversed. Therefore the optimal current protocol is such that it 
exactly time-reverses the purely relaxational trajectory of an isolated system. Such a 
relaxation trajectory must be thought of as starting at the energy maximum (i.e. 
the switching point) and winding down towards the energy minimum at i? = 0. This 
is a rather general statement which is based only on the assumption that the energy 
loss is quadratic in the spin-current (the latter is very nearly true practically in all 
systems investigated so far). 

To put this observation on a more practical level we notice that the energy- 
dependent critical current (|1.39p appears to be roughly energy independent. Equation 
(jl.42p implies then that the optimal current is very nearly a constant given by twice 
the critical switching current Ic- This is shown in Fig. ll.Gl via Monte Carlo simulations 
of the Joule heating as a function of the spin current for various confidence levels. At 
low confidence levels < 2Ic due to large fluctuations, however at confidence levels 
nearing 100% we see the predicted I°p* « 2Xc. 

While we have established that the optimal current pulse is approximately rect- 
angular with the amplitude sa 2Ic, we have not discussed yet its duration. 
As discussed above the distribution of the switching times is well described by the 
equilibrium thermal distribution of initial energy oc e~^""/'^^^ and the determinis- 
tic duration ()1.40|) . Thus the times needed to reach the switching point, which are 
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Fig. 1.6 Simulated Joule iieating during the switciiing process in arbitrary units vs. spin 
current for various confidence levels p. From top to bottom p = 0.99, p = 0.5, p = 0.2. Here 
a = 0.01, = 0.028Ms, = and Eo = WOfcsT. 

oc a^^ ln(i?o/i?ini), where Eq = iiqH^Ms/2 is the hight of tlie energy barrier, are also 
scattered accordingly. As a result, any pulse of a finite duration achieves switching 
only with a certain probability p < 1. If this probability, i.e. error tolerance, is speci- 
fied such that 1 — p ^ 1, all realizations with ii^ini > kBT[l — p) should undergo the 
switch. This dictates that the duration of the optimal current pulse scales as 



opt 



In 



En 



a \ kBT{l — p) 



(1.44) 



where to — (jMs) ^ . This result is shown in Fig. Il.Tl for higher switching probabilities 
and is in a good agreement with the simulated data. 




Fig. 1.7 Simulated optimal switciiing time as a function of tlie confidence level (crosses), 
eqn ifLii)! (full line). Here a = 0.01, = 0.028Afs, Hj^ = Afs, Eo = UOksT and Is = XJ"'. 



One can imagine sufficiently complex switching protocols, utilizing the preccssional 



Spectral width of the steady state precession 17 



frequency to achieve an optimal path. For example, Ref. ( |Cui et ai, 2008[ ) used current 
pulses resonant with the orbit period to excite the system. In the realm of relatively 
low frequency protocols, where the averaging over the fast angle is justified, twice the 
critical current pulse provides the best approximation to the optimal strategy. The 
duration of this pulse weakly depends on the confidence level and can be estimated 
according to cqn (|1.44|) . 

1.6 Spectral width of the steady state precession 

Another important application of spin-torque devices is provided by spin-torque oscil- 
lators. These devices arc similar to spin valves, discussed above, but rather rely on the 



steady state precession to generate microwave signal ( Deac et ai, 2008 Mizukami et ai, 2001 



Krivorotov et al., 2005||Nazarov et al., 2006|[Pribiag et ai, 2009[ [Georges et al, 2009 



Ruotolo et ai, 2009| [Cheng et ai, 2010| [Finocchio et ai, 2010[ [Rippard et ai, 2010[ 



Finocchio et al, 2011 1. Because of the high precision needed in these applications it is 



important to understand how the noise can impact spectral line shape of the oscillator. 
There is a number of experimental measurements of the irradiated frequency spectrum, 
which indicate a non-monotonous dependence of the spectral width on various pa- 
rameters such as the spin-current amphtude ( Kiselev et al., 2003[ [Mistral et ai, 2006[ 
[Tiberkevich et ai, 2007[ [Georges et al, 2009p . 

The steady state magnetization precession (SSMP) occurs along SW orbit when the 
energy pumped by the spin-torque is compensated by the dissipation. If this is the case, 
the time derivative of the total energy along a certain orbit, eqn (jl.l3p . vanishes, E = 
dut E = E, full line in Fig. 11.41 Generically, SSMP is realized within a certain range of 
the spin-currents below the critical one Is ^ Ic- The corresponding energy E = E{Is) 
of the stable SW orbit typically runs from zero to the switching energy Eq, with the 
increase of the spin-current. The energy distribution corresponding to the steady state 
magnetization precession can be obtained from eqn p.38p . As follows from eqn p.l3|l 
the generalized force E = — -Fg -I- Is ■ = 0, thus the energy distribution can be 



obtained by expanding the force as i? ~ {E — E). Then, employing cqn (|1.38p 

and neglecting the energy dependence of the diffusion coefficient, one obtains 



^SSMP 




(E) (X exp 



2Vt 



(1.45) 



Therefore, the SSMP energy distribution is approximately Gaussian, centered at the 
energy E with the width determined by the diffusion coefficient and the derivative of 
the generalized force E over the energy. 

To investigate the spectral width of the precession spectrum, one needs to restore 
the stochastic equation (|1.17p for the angular variable (p{t) running along SW or- 
bits. Employing the fact that the correlator (|1.3ip of the stochastic magnetic fields is 
isotropic, it may be written as 

I^E I [iVl X JrloffJ I 

where ^^{t) is a scalar Langevin force with the correlator {^ip{t)£,ip{t')) = 2D(M.)6{t — 
t'). It is thus the function \g^p{E,ip)\, which determines the frequency spectrum. Some 



18 Dynamics of nano-magnetic oscillators 



of it universal features can be inferred without relying on a specific model. In particular, 
it can be concluded that the function \g^{E,ip)\ diverges as E tends to zero. Indeed, 
at small energies E, the characteristic size of SW orbit goes to zero (at zero energy 
the orbit shrinks to a point), while the influence of noise on the change of angle ip 
grows. The character of the singularity can be obtained by considering the simplest 
model of steady state precession in a system with an easy axis anisotropy and both 
external magnetic field and spin current being parallel to the easy axis. In this case the 
energy is a function of cos 9 and SW orbits are circles, while ip is the azimuthal angle 
and \gtp{E,(p)\ cx l/sin0 ( [Chudnovskiy et ai, 2008 1. Taking into account that the 



energy minimum corresponds to = 0, and expanding at small angles E{cos9) cx 
one concludes that |g^(i?)| cx 1/Ve. Therefore, at small energies of SW orbits, the 
noise strength decreases with the growth of the energy. This feature is responsible 
for the decreasing linewidth of the microwave spectral power with the increase of the 



spin-current (and thus increase of E), as discussed in Refs. (Tiberkevich et ai, 2007 



|Kim et ai, 2008D assuming the equilibrium thermal noise ( |Brown, 1963[ ). 



To evaluate the spectral power, one notices that the precessing magnetization 
induces the oscillating magnetic field H in a waveguide, which is a periodic func- 
tion of the angle cp, i.e 'H(<p) = 'H(<p + 27r). The spectral power is given by S{u}) = 
J dt {'H{(p{t))'H{(p{0))) ^ g-jwt g^^-^^ according to eqn (|1.46p it consists of series of peaks 
at multiple integers of J^^. Focusing on the linewidth of the n-th harmonics oj = n^}j^, 
one needs to calculate 



5(w)cx y di(e™[^(*)-'^(")l)^ e-'^"*; cp{t) ^ Qst + J \g^{E,ip)\^^it')dt' , 

(1.47) 

where we have employed the formal solution of eqn (|1.46p . The averaging over the 
noise is thus reduced to the Gaussian integral, which gives 

S{lu) <x J dt cxp l^tinQE - t^)t - y ^ D(M{t')) \g^{E, <f{t'))\^dt'^ . (1.48) 

This expression should be now averaged over the stationary energy distribution (|1.45p . 
Below we shall assume that the latter is rather narrow and forces the energy to be 
close to the stationary value E. In this case the spectrum is Lorentzian centered around 
u! = nftj^ with the width given by n^A^, where is the zeroth Fourier harmonics 
of the last term in the exponent in eqn (jl.48p 

^^Jo ^eJo |[MxHcff]| 

Finally, using the equation of motion (|l.ip . one may rewrite the spectral width as 
an integral along the static SW orbit with the energy E 

A.-^/^(M) MMxH.,].M 

iPl J ^ ' |[MxHeff']|' 

compare this result with that for the energy diffusion coefficient ()1.24|) . Proportion- 
ality of the spectral width to the noise correllator £)(M), eqn (|1.3ip . implies that 
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Voltage V 



Fig. 1.8 Calculated linewidth of the microwave spectral power vs. applied voltage in units of 
temperature. Inset: calculated microwave power (solid line) and the average precession angle 
(dashed line) vs. voltage. Parameters: Ms/{h'y) = 10, K'yHl/kB = 3K, T = IK, ao = 0.01. 
Conductances in units /h: Gp = 0.181, Gap = 0.019. dlsi/dV = O.Ole. 

it is proportional to temperature if T > eV . On the other hand, at small temper- 
ature and/or larger voltage V the correlator is dominated by the non-equilibrium 
noise and one expects A cx 1^ cx I^- Therefore one expects decrease of the spec- 
tral width with decreasing temperature and saturation on a voltage-dependent level 
dSankey et al, 2005D . 

Figure rOl shows the calculated spectral linewidth as a function of the applied volt- 
age at T — IK. Notice the non-monotonous behavior, similar to the one observed in a 
number of experiments ( Rippard et ai, 2004[|Sankey et al., 2005 Rippard et al., 2006 



Mistral et al, 2006 Tiberkevich et al., 20071. The initial decrease of the spectral width 



is associated with the angular factor jM x Hofi|^ (x E in the denominator of eqn (|1.49p 
and the fact that the stable energy E grows with the applied spin-current. On the other 
hand, the subsequent reversal of the trend and growth of the spectral width is due 
to the fact that £>(M) grows with the voltage, because of the spin-torque shot noise, 
eqn (fOTj) . 



1.7 Conclusion and acknowledgments 

Miniaturization of spintronics devices puts the stochastic effects in the magnetization 
dynamics on the cutting edge of investigations. While the fundamentals of stochastic 
description of magnetization in terms of the FP equation have been developed decades 
ago, the highly nonequilibrium regimes used in modern devices require a specific de- 
scription of non-equilibrium noise. In this work, a step towards such a description has 
been reviewed. As we showed, despite a large variety of dynamical regimes, they all 
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have in common the time scale separation between the fast and almost energy con- 
serving magnetization precession and the slow energy change. This point allows to 
develop a theoretical description of magnetization dynamics in terms of slowly vary- 
ing energy of the precessional orbit. The resulting FP equation is one dimensional, 
which greatly simplifies its treatment in comparison to the initial FP equation for the 
magnetization vector. In this work we applied such an approach to the analysis of 
current induced magnetization switching and steady state magnetization precession, 
obtaining results on the optimization of the current switching protocol and on the 
linewidth of radiation emitted by steady state magnetization precession. Further de- 
velopment of the energetic description is required by the dynamical regimes, where 
resonances between the precessional orbits with different energies are induced by ap- 
plying an ac- modulated spin current as reported recently in e.g. ( [Cheng et al., 2010[ ), 
as well as for the description of spin switching by spin current pulses of complex form 
( Cui et al, 20081 ). 

We are grateful to Ilya Krivorotov for fruitful discussions. ALC acknowledges sup- 
port from DFG through the Collaborative Research Center 668. AK and TD were 
supported by NSF Grant DMR-0804266 and U.S.-Israel Binational Science Founda- 
tion Grant 2008075. 
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